Reduced carbon solubility in Fe nano-clusters and implications for the growth of 

single-walled carbon nanotubes 
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Various diameters of alumina-supported Fe catalysts are used to grow single-walled carbon nan- 
otubes (SWCNTs) with chemical vapor decomposition. We find that the reduction of the catalyst 
size requires an increase of the minimum temperature necessary for the growth. We address this 
phenomenon in terms of solubility of C in Fe nanoclusters and, by using first principles calculations, 
we devise a simple model to predict the behavior of the phases competing for stability in Fe-C nano- 
clusters at low temperature. We show that, as a function particles size, there are three scenarios 
compatible with steady state-, limited- and no-growth of SWCNTs, corresponding to unaffected, 
reduced and no solubility of C in the particles. The result raises previously unknown concerns about 
the growth feasibility of small and very-long SWCNTs within the current Fe CVD technology, and 
suggests new strategies in the search of better catalysts. 
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Among the established methods for single-walled car- 
bon nanotubes (SWCNTs) synthesis [J 0, the low 
temperature catalytic chemical vapor decomposition 
(CCVD) technique is more appropriate for growing nan- 
otubes on a substrate at a target position and, thereby, 
can accelerate the integration of this unique material in 
hybrid electronics. Reported synthesis of SWCNTs by 
CCVD at reactor temperature as low as ~350-450°C was 
achieved by using hydrocarbons with extremely exother- 
mic catalytic decomposition reaction [J 0], which may 
release significant free energy (e.g. for acetylene 260 
kJ/mol at 650° C @), affecting the temperature of the 
catalyst under certain conditions. 

Considering the vapor-liquid-solid model (VLS) as the 
most probable mechanism for S WCNT growth [1, [13, EH , 
an alternate approach for lowering the growth temper- 
ature can be the reduction of the catalyst size, since 
the Gibbs-Thompson model predicts a decrease of the 
melting temperature with decreasing cluster size 15l. Il6, 
17j and the synthesis temperature is correlated to the 
catalyst-carbon melting and eutectic points [l8, 1^, 20] ]. 
While small catalyst particles nucleate small diameter 
tubes (this may lead to the diameter control growth, an- 
other fundamental problem that hinders nanotube ap- 
plication) , they also affect the morphology of the formed 
carbon structures 21 1 , the kinetics of the growth 22 , 2^ , 
and the solubility of carbon available for the growth pro- 
cess, which requires understanding of the widely debated 
question about the property of the catalyst particles, as- 
sociated with the fundamental problem of their thermo- 
dynamic state 151 ] . The main challenges are the un- 



exposed peculiar features of small catalysts and their 
binary phase diagrams with carbon. In this work, we 
address the carbon solubility problem in Fe nanoparti- 
cles by studying size-dependent growth of SWCNTs by 
CCVD, and the thermodynamics of competing phases in 
Fe-C clusters with ab initio modeling. 



In our experiments, Fe catalysts supported on alumina 
powder were prepared using a the common impregnation 
14j |. To find the minimum synthesis tempera- 
as a function of catalyst size, different par- 



method 
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tides' dimensions were obtained by varying the Fe:Al20 



24j corresponding 
1.4±0.7, 



molar ratio 1:15, 1:25, 1:50 and 1:100, 
to particles of diameters ~ 3 ± 0.6, ~ 2 ± 0.8, 
and ~ 1.3 ± 0.7 nm, respectively. The growth of SWC- 
NTs was performed with CCVD at temperatures between 
650-900°C, for 90 min over the pre-reduced catalysts [25j |. 
To avoid increasing the catalyst temperature, methane 
was used as the carbon source since the formation of 
SWCNTs from methane is not exothermic. For each case 
T^Jynth was determined by analyzing the Raman spectra 
(Thermo Nicolet Almega Raman spectrometer equipped 
with a CCD detector, laser excitations of 532 and 785 
nm) of the samples collected after synthesis at various 
temperatures. 
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FIG. 1: Evolution of minimum synthesis temperature 
(T^th) for growing SWCNTs with Fe:Al 2 3 molar ratio. 
The inset shows the Raman spectra for samples obtained with 
a catalyst molar ratio of 1:50 at 820, 850 and 865°C. 
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Figure Q] shows the T™" t/l dependence on the Fe:Al 2 3 
molar ratio. The inset presents a typical example of Ra- 
man spectra after the synthesis between 820 and 865°C. 
As can be seen from molar ratio of 1:50, SWCNTs can 
be grown at 865° C but not at 850° C even though car- 
bon deposition is observed. In this case, T™"^ was es- 
timated to be -865±9°C. We find that T 



synth 



increases 



with decreasing catalyst size (molar ratio), contrary to 
what may be expected from the Gibbs-Thompson model 
(lij ]. This observation indicates that decomposition of 
the hydrocarbon alone is not enough to grow nanotubes 
and that the temperature must be increased to ensure 
that a certain amount of carbon dissolves into the parti- 
cle (considering that the maximum solubility of C in Fe 
depends on the catalyst size, as shown later). In fact, 
temperature must be increased to overcome the loss of 
solubility of C in the catalytically-active phase compet- 
ing for thermodynamical stability with a nucleating car- 
bide, and not only to enhance diffusion of C (otherwise 
below T^nth we would have shorter nanotubes, instead 
of their absence). 

We believe that the origin of this apparent paradox lies 
in a novel phenomenon, i.e., a reduced solubility of C in 
Fe nanoparticles. Within the VLS framework with bulk 
diffusion as the rate- limiting step [ll|, 0, O EH? this 
implies an increase of temperature to achieve compara- 
ble amount of dissolved carbon to allow growth. In Ref. 
[lH we have shown that the eutectic point (x^ ut ,T euf ) of 
Fe-C clusters shifts toward lower carbon concentrations, 
x eut> with decreasing particle size (Figure 8 of [lit]). Due 
to the high energetic cost for bringing bulk cementite off- 
stoichiometry (in the Fe-C phase diagrams FesC forms 
two-phase regions with austenite (7) and ferrite (a) with- 
out going off-stoichiometry (51|), the most probable cause 
of the shift of x^ ut is a reduced solubility of C. This does 
not necessarily imply that the total amount of C present 
in the nanoparticle decreases, because solubility counts 
only the carbon dissolved randomly in the solid Fe-rich 
phases, (in equilibrium with FeaC (below T eut ) or with 
the liquid/ viscous phase (above T eut )). 

The accurate analysis of the phenomenon can be 
achieved by calculating the interplay between the phases 
competing for stability at the temperatures of the pro- 
cess. For nanoparticles, this task is generally unsolvable, 
although qualitative information can be extracted from 
approximate zero-temperature first-principles modeling. 
In such approaches, by comparing the formation energies 
of the candidate phases, we determine the stability of 
the system at low-T and give indications for higher tem- 
perature behavior. The ab initio simulations presented 
here are performed with VASP [28j |. using projector aug- 
mented waves (PAW) [29j and exchange-correlation func- 
tionals as parameterized by Perdew, Burke, and Ernzer- 
hof (PBE) 30] for the generalized gradient approxima- 
tion (GGA). Simulations are carried out with spin polar- 
ization, at zero temperature, and without zero-point mo- 



tion. All structures are fully relaxed. Numerical conver- 
gence to within about 2 meV/atom is ensured by enforc- 
ing a high energy cut-off (500 eV) and dense k-meshes. 
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FIG. 2: (color online). Size-pressure approximation for Fe 
nanoparticles. Given a spherical particle, we estimate the hy- 
drostatic pressure due to the surface curvature by calculating 
the deviation of the average bond length inside the cluster, 
Ad nn = d% n — d nn , and comparing it to the deviation of the 
bond length for the bulk material as a function of hydrostatic 
pressure. 



Our model is based on the following assumptions: (i) 
Mechanism. The behavior of carbon is determined by the 
interplay of four competing phases as a function of cat- 
alyst size: pure bcc-Fe, C dissolved in ferrite (a-FeCV), 



ordered cementite (FeaC), and carbon SWNTs [14, 
The pure-Fe phase is taken to be bec because our simu- 
lations are aimed to explore the low temperature regime 
of catalytic growth. The a-FeC x phase is simulated by 
taking samples of bec supercells with different concen- 
trations of interstitial carbon (Fe32C, Fe24C, Fe^C). We 
are not required to generate truly random phases by using 
the special quasirandom structure formalism, because the 
low concentration of carbon, guarantees enough distance 

" 32|. In 



between replica of C to achieve convergence [31 . 
addition, higher concentrations of C are not required to 
be explored, because even in bulk a-¥eC x the solubility 
is small [2rj | . (ii) Carbon source. Free carbon atoms come 
from the dissociation of the feed-stock on the surface of 
pure-Fe and random FeC^ catalysts only. Formation of 
cementite stops the process due to its different activity 
and diffusion properties (as show in Fig. 1 of 14j and ref- 
erences therein), (iii) Nanotubes diameter. To minimize 
the curvature energy of the tube, active catalysts produce 
nanotubes that have a similar diameter as the particle 
(CVD experiments have shown clear correlation between 
the two diameters 27]). (iv) Size-pressure approxima- 
tion. In nanoparticles, surface curvature and superficial 
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dangling bonds are responsible for internal stress fields 
which modify the atomic bond lengths inside the parti- 
cles. For spherical clusters, the induced stress fields are 
hydrostatic and the compressions are isotropic. The phe- 
nomenon can be modeled with the Young-Laplace equa- 
tion Ap = 2j / R where the proportionality constant 7 
(surface tension for liquid particles) can be calculated 
with ab initio methods. As a first approximation, by 
neglecting all the surface effects not included in the cur- 
vature, the study of phase diagram for spherical parti- 
cles can be mapped onto the study of phase diagrams for 
bulk systems under the same hydrostatic pressure pro- 
duced by the curvature, as depicted in Fig. [2] 33]. It is 
worth mentioning that our reference for carbon is taken 
to be the zero pressure nanotube phase, different from 
the other carbon references used for investigating Fe-C 
under pressure 34 , 35|, 36| . 

Figure [2] shows the implementation of the "size- 
pressure approximation" for Fe nanoparticles. On the 
left hand side we show the ab initio calculations of the 
deviation of the average bond length inside the cluster 
Ad nn = d„ n — d nn (d„ n = 0.2455 nm is our bulk bond 
length), for bcc particles of size N = 59, 113, 137, 169, 307 
and 00 (bulk) as a function of the inverse radius (l/R)- 
The particles were created by intersecting a bcc lattice 
with different size spheres. The particle radius is de- 
fined as l/R = 1 /N scp j 1 /Ri where the sum is taken 
over the atoms belonging to the surface convex polytope 
(N scp vertices) and FL are the distances to the geometric 
center of the cluster |38[ . The left straight line is a linear 
interpolation between 1 / R and Ad nn calculated with the 
constraint of passing through l/R = and Ad nn = 
(N = 00, bulk). The right hand side shows the ab initio 
value of d nn in bulk bcc Fe as a function of hydrostatic 
pressure, P. The straight line is a linear interpolation 
between P and Ad nn calculated with the constraint of 
passing through P = and Ad nn = (bulk lattice). 
By following the colored dashed paths indicated by the 
arrows we can map the analysis of nanoparticles' stabil- 
ity as function of R onto bulk stability as function of 
P, and obtain the relation between the radius of par- 
ticle/nanotube and the effective pressure PR— 2.46 
GPa-nm. It is important to mention that our 7 = 1.23 
J/m is not a real surface tension but an ab initio fit- 
ting parameter describing size-induced stress in nanopar- 
ticles. In addition, it compares well with the experimen- 
tal value of the surface tension of bulk Fe at its melting 
point ~ 1.85 J/m 2 [37|. 

Figure [3] shows the evolution of competing phases as 
function of concentration and particle's size. In panels (a) 
through (d), we plot formation energies, Ef [Fei- y C y ] = 
£[Fei_j,Cj,] - (1 - y)E[Fe] - yE[C], calculated with re- 
spect to the pure constituents of the reaction, bcc-Fe and 
carbon-SWCNTs (with the same diameter as the parti- 
cle). The dashed green lines interpolating FeigC through 
Fe, are used to estimate the formation energies of the 
random phase with maximum solubility of carbon (a- 
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FIG. 3: (color online). The scenarios described in the text 
with their implication in CVD SWCNT growth with Fe nano- 
catalysts. P is the pressure at which the calculations are per- 
formed, and R is the radius of the corresponding SWCNT 
within the size-pressure approximation. E/[Fe3C]=E/[a- 
FeCttooicw] for R mm ~ 0.58 nm (P ~ 4.3GPa). 

FeC^g qq 102 , black squares near the origin in Fig. [3] la- 
beled as FeC*) corresponding to 0.022 weight% [26^ A 
structure at a given composition is considered stable (at 
zero temperature and without zero-point motion) if it 
has the lowest formation energy for any structure at this 
composition, and, if on the binary phase diagram, it lies 
below the convex hull of tie lines connecting all the other 
stable structures [39j, |4fJ • Phases lying above the convex 
hull and with small positive formation energies, might be 
explored by the thermodynamics of the system through 
configurational and vibrational entropic promotion. 

By varying the radius of the particle, the stability of 
the competing phases, a-FeC^ and Fe 3 C, changes con- 
siderably. There are three possible scenarios. Scenario 
I is shown in panels (a) and (b). For big particles, 
R» R m in, Fe$C has formation energy higher than the 
max solubility phase (a-FeC^g 00102)- Therefore the pol- 
lution of carbon at low and medium temperature, can- 
not cause the big particle to undergo phase transition 
by nucleating cementite. Hence, such particles remain in 
the catalytically-active random oFeC x state, by kee ping 
their concentration of carbon between to ~ 0.102% 15 1 
(solubility is unaffected), and by implying a balance be- 
tween in- and out- flows of carbon which can guarantee 
the steady state growth of nanotubes. Thermodynam- 
ically, in this regime, SWCNTs, MWCNTs and carbon 
fibers could be grown indefinitely and the only limitation 
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is the availability of carbons feed-stock 0. Thus, ex- 
periments performed with particles of these sizes would 
be described by Arrhenius equations governing catalytic 
activity and diffusion properties. The minimum radius, 
Rmin, can be estimated by interpolating the pressure at 
which the energy of FesC and a-FeC^g. 00102 are equal, 
and mapping such pressure in the size-pressure relation 
of Fig. H We obtain P ~ 4.3GPa and R mm ~ 0.58nm. 
Scenario II is shown in panel (c). For particles of size 
R~R 

mim the max solubility phase and Fe3C have simi- 
lar formation energies. The energetically competing FeaC 
causes depletion of C in a-FeC (reduced solubility), nu- 
cleation of ordered FesC, and overall reduction of catalyt- 
ically active random Fe. If exposed to hydrocarbons at 
elevated temperatures, such particles would be capable 
of dissociating carbon and growing SWCNTs with con- 
comitant nucleation of cementite. Such nucleation slowly 
poisons the particle and terminates the growth. In this 
regime, SWCNTs can be produced up to a certain critical 
length depending on the net flow of carbon. By varying 
R, the critical length goes from infinity (R > R m i n ) to 
zero (R < R m in)- Experiments performed with particles 
of size R ~ R m in would show on/off growth at low tem- 
perature and Arrhenius behavior at high temperature. 
Scenario III is shown in panel (d) . By further reducing 
the size of the particle, R << Rmin, the formation en- 
ergy of cementite becomes negative (Ef{a-FeC^Q .00102] 
is always very close to zero). The stability of FesC over 
the range to 25% carbon in the phase diagram, indi- 
cates that the nucleation of FesC occurs simultaneously 
with the carbon pollution. The particle transforms into 
FesC as rapidly as the availability of feed-stock allows, 
thus no random phase coexists at low temperature (max 
solubility of C is zero), and no out-flow of carbon occurs. 
Particles with R < R m i n cannot grow SWCNTs, and 
Rmin can be considered as a lower limit for SWCNTs' size 
in low-temperature CVD growth with Fe nano-catalysts. 
Experiments performed with such particles would result 
in FesC nanoparticles and no appreciable nanotube pro- 
ductions. 

We address the extension of the model to higher tem- 
peratures in a qualitative framework 41(. Increasing 



the temperature of the reaction and considering fee Fe 
nanoparticles has two consequences. FesC becomes sta- 
ble at medium temperature [26j and, by reducing particle 
size, the stability versus size follows similar arguments as 
the case of the bec reference. Thus R m in — > 00 and a 
steady state SWCNT growth is, a priori, unobtainable. 
However, since the maximum solubility of C in austenite 
7-FeCa; is bigger than in ferrite a-FeC x (due to the fact 
that interstitial holes in the fee lattice are larger than 
those in the bec lattices (~2.0 wt% in j-FeC x [26:])), the 
life of the catalytic particle can somehow be longer than 
the previous scenarios 2 and 3, and the grown SWCNTs 
can be quite long. The lack of steady state at high tem- 
perature can be addressed by alloying the nanoparticle 
with other metals to reduce cementite stability and by 



simultaneously promoting other more stable alloyed ran- 
dom phases with considerable catalytic activity. 

In this manuscript we have addressed the thermo- 
dynamic limits of bec Fe nano-catalysis technology for 
SWCNT growth with low temperature CVD. Using a 
"size-pressure approximation" and zero temperature ab 
initio modeling we identify ranges of nanoparticle sizes 
which are compatible for steady state-, limited- and no- 
growth of SWCNTs corresponding to unaffected, reduced 
and no solubility of C in the a-FeC^, nanoparticles. 
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